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ABSTRACT 

We have obtained radial gradients in the spectral features A4000-A break (D4000) and 
Mg2 for a sample of 11 central cluster galaxies (CCGs): 8 in clusters with cooling flows 
and 3 in clusters without. After careful removal of the emission lines found within the 
C*~) ' D4000 an d the Mg2 bandpasses for some objects, the new data strongly confirm the 

\ correlations between line-strength indices and the cooling flow phenomenon found in 

our earlier study (Cardiel, Gorgas & Aragon-Salamanca 1995). We find that such cor- 
relations depend on the presence and characteristics of emission lines in the inner 
regions of the CCGs. The nuclear indices are correlated with the mass deposition rate 
(M) only when emission lines are found in the central regions of the galaxies. The 
\ central D4000 and Mg2 indices in cooling flow galaxies without emission lines are com- 

■^j- ■ pletely consistent with the indices measured in CCGs in clusters without cooling flows. 

f^S CCGs in cooling flow clusters exhibit a clear sequence in the D4ooo _ Mg2 plane, with 

OO . a neat segregation depending on emission-line types (Heckman et al. 1989) and blue 

0^ ' morphology (McNamara 1997). This sequence can be modelled, using stellar popula- 

tion models with a normal IMF, by a recent (~ 0.1 Gyr old) burst of star formation, 
although model uncertainties do not allow us to completely discard a continuous star 
formation or a series of bursts over the last few Gyr. In CCGs with emission lines, the 
gradients in the spectral indices are flat or positive inside the emission-line regions, 
suggesting the presence of young stars. Outside the emission-line regions, and in cool- 
^ ' ing flow galaxies without emission lines, gradients are negative and consistent with 

those measured in CCGs in clusters without cooling flows and giant elliptical galaxies. 
Index gradients measured exclusively in the emission-line region correlate with M . 
k>( , Using the same population models we have estimated the radial profiles of the mass 

. transformed into new stars. The derived profiles are remarkably parallel to the ex- 

pected radial behaviour of the mass deposition rate derived from X-ray observations. 
Moreover, a large fraction {probably most) of the cooling flow gas accreted into the 
emission-line region is converted into stars. In the light of this new data, we discuss 
the evolutionary sequence suggested by McNamara (1997), in which radio triggered 
star formation bursts take place several times during the lifetime of the cooling flow. 
We conclude that this scenario is consistent with the available observations. 

Key words: stars: formation - cooling flows - galaxies: elliptical and lenticular, cD 
- galaxies: starburst - galaxies: stellar content - X-rays: galaxies 



1 INTRODUCTION 

The cooling time of the hot intra cluster medium is shorter 
than the Hubble time in the central regions of most 
galaxy clusters (e.g. Arnaud 1988, Edge, Stewart & Fabian 
1992), which suggests that large amounts of cool material 
should be settling at the cluster centre. This process, nor- 
mally known as a cooling flow, may take place at rates 
of M ~ 100 M yr" 1 (Fabian 1994), which implies that 



A^Totai ~ 10 Mq of gas would cool and accrete onto the 
central cluster elliptical galaxies over the lifetime of the clus- 
ters. Considering the idealised picture of an homogeneous 
cooling flow, all the material will be deposited in the cen- 
tre of the cluster. However, X-ray data show that the gas is 
highly inhomogeneous and, very likely, most of the gas must 
be cooling out at large radii, whilst only a fraction reaches 
the centre of the flow. In fact, mass deposition profiles are 
reasonably fitted by M(< r) cx r (Thomas, Fabian & Nulsen 
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1987), although some clusters do not follow this trend be- 
yond a certain radius r ^ 60 kpc (David, Jones & Forman 
1994; Irwin & Sarazin 1995). For comprehensive reviews on 
cooling flows see Fabian, Nulsen & Canizares (1984, 1991), 
Sarazin (1986) and Fabian (1994). 

A question of open debate is the nature of the final 
repository of the cool gas. The discovery of significant excess 
absorption in the soft X-ray spectra of some galaxy clusters 
containing cooling flows (White et al. 1991; Johnstone et al. 
1992; Miyaji et al. 1993; Allen et al. 1993, 1995; Fabian et al. 
1994; Irwin & Sarazin 1995; Allen & Fabian 1997a, 1997b) 
requires ~ 10 11 -10 12 M© of absorbing cold gas, probably in 
small, pressure-confined, magnetised clouds (Daines, Fabian 
& Thomas 1994). This value is comparable to A/xotai, the 
mass deposition rate integrated over the age of the clus- 
ters. Fabian, Johnstone & Daines (1994) have postulated 
that dust can form in the cold clouds embedded in cooling 
flows, which would explain why the clouds are almost unde- 
tectable outside the X-ray wavelength range. However, there 
is a great uncertainty concerning the nature of the cold ma- 
terial producing this absorption (Sarazin 1997, Laor 1997). 
Voit & Donahue (1995), from the limits imposed by 21 cm 
and CO observations (McNamara, Bregman & O'Connell 
1990; Antonucci & Barvainis 1994; McNamara & Jaffe 1994; 
O'Dea et al. 1994; O'Dea, Gallimore & Baum 1995), suggest 
that it is unlikely that significant amounts of cold gas can 
remain undetectable. On the other hand, very recently Jaffe 
& Bremer (1997), using 7f-band spectroscopy, have found 
that the inner few kpc of central cluster galaxies (CCGs) in 
cooling flows exhibit strong emission in the H2(1-0)S(1) line, 
which is not seen in a comparable sample of non cooling flow 
galaxies. These authors suggest that it is likely that a large 
mass of additional molecular material covers the cooling flow 
as a whole, producing the soft X-ray absorption. 

Since it seems clear that at least a fraction of the 
cool gas must reach the inner regions of the CCGs, it is 
reasonable to think that part of this gas could accumu- 
late as molecular clouds which, presumably, would form 
stars. Unusual diffuse blue light and colours (Sarazin & 
O'Connell 1983; Bertola et al. 1986; Johnstone, Fabian & 
Nulsen 1987 —hereafter JFN87— ; Romanishin 1987; Craw- 
ford et al. 1989; McNamara & O'Connell 1989, 1992 - 
hereafter M089, M092— , 1993; Thuan & Puschell 1989; 
N0rgaard-Nielsen, J0rgensen & Hansen 1990; Allen et al. 
1992; Crawford & Fabian 1993; Crawford et al. 1995; Allen 
1995; McNamara 1995; Hansen, J0rgensen & N0rgaard- 
Nielsen 1995; Cardiel, Gorgas & Aragon-Salamanca 1995 
— hereafter CGA95 — ; Melnik, Gopal-Krishna & Terlevich 
1997), strong low- ionization optical line emission (Baum 
1992, Donahue & Voit 1997, and references therein), and 
powerful radio sources (see Jaffe 1992, Burns et al. 1997, 
for reviews) are usually associated with the inner few kpc of 
central dominant galaxies in clusters with clear X-ray sig- 
natures of cooling flows. Although with a large scatter, the 
strength of the optical anomalies are correlated with mass 
deposition rates (JFN87; M089; M092; Heckman et al. 1989 
—hereafter HBvBM— ; CGA95; McNamara 1997). In most 
cases star formation has been invoked to reproduce the ex- 
cess blue light in the continuum, with typical derived star 
formation rates ranging from a few to some tens of solar 
masses per year. It is important to stress that these values 



normally only account for ^ 10 per cent of the total X-ray 
derived mass deposition rates. 

McNamara & O'Connell (1993) found that the blue 
colours observed in A 1795 and A 2597 are located in patches 
that are coincident with the radio lobes of the central galax- 
ies, which suggests a connection between the radio source 
and the star formation (or, at least, the source of the blue 
light). Although power-law spectra could also be an alter- 
native explanation for the origin of the excess blue light 
(e.g. Crawford & Fabian 1993), the absence of polarization 
in (7-band observations of the blue lobes in A 1795 (Mc- 
Namara et al. 1996b) makes this scenario unlikely. In ad- 
dition, observations of this object with the Hubble Space 
Telescope (McNamara et al. 1996a; Pinkney et al. 1996), 
and the Ultraviolet Imaging Telescope (Smith et al. 1997), 
strongly support that the origin of the blue lobes could be 
an accreted population formed from the cluster cooling flow 
or from episodes of star formation induced by accretion of 
stripped material from several gaseous dwarf galaxies. 

In this work we address the possible connection between 
star formation and the cooling flow phenomenon through the 
study of the stellar populations in CCGs. For this purpose 
we analyse line-strength gradients in Mg2 (Faber et al. 1985) 
and the A4000-A break (D4000, definition in Bruzual 1983) 
in a sample of 11 brightest cluster galaxies. This paper is an 
extension of a previous work presented in CGA95 after the 
inclusion of extensive new spectroscopic data. 

The galaxy sample, reduction of the data and er- 
ror analysis are presented in Sections 2 and 3. The nu- 
clear line-strength and gradients are described in Sec- 
tion 4. The discussion of these results is presented in Sec- 
tion 5, and our conclusions are summarised in Section 6. 
Ho = 50 km s^ 1 Mpc -1 has been assumed throughout the 
paper. 



2 GALAXY SAMPLE AND OBSERVATIONS 

We have obtained long-slit spectra of CCGs in two runs with 
different telescopes. Since the aim of the observations was 
to get reliable line-strength gradients, they were performed 
on dark nights. The first run (August 1994) was carried 
out with the 3.5-m telescope at the German- Spanish As- 
tronomical Observatory at Calar Alto (Almeria, Spain). We 
used the Cassegrain Twin Spectrograph with a blue coated 
TEK CCD in the red channel. In the second run (December 
1995) we observed with the 4.2-m WHT at the Roque de los 
Muchachos Observatory (La Palma, Spain), using the blue 
arm of ISIS with a blue coated TEK CCD. A summary of the 
observations, including relevant instrumental parameters, is 
given in Table 1. 

The sample of CCGs was chosen to cover a wide range 
of mass deposition rates and complements our previous data 
(CGA95). It includes 8 galaxies in clusters with cooling flows 
and 3 galaxies in clusters without cooling flows. A descrip- 
tion of the sample, including exposure times and position 
angles of the spectrograph slit, is listed in Table 2. We have 
re-observed the central galaxy of A 644 since previous mea- 
surements (CGA95) were not clearly understood. We also 
observed G-K giant stars to guarantee that our Mg2 indices 
are in the appropriate photometric system (see below). 
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3 DATA REDUCTION AND ERROR 
ANALYSIS 

Reduction of the data was carried out using our own reduc- 
tion package REDWrE (Cardiel & Gorgas 19980). We fol- 
lowed a standard reduction for long-slit spectroscopic data: 
bias and dark subtraction, cosmic ray cleaning, flat-fielding, 
wavelength calibration, C-distortion correction, sky subtrac- 
tion, extinction corrections, centring and binning of the 
spectra. The process is similar to that explained in CGA95, 
and its description is not repeated here. However, we present 
a summary of the error sources and the procedure followed 
for their estimation. 



3.1 Random errors 

We describe below the main sources of random errors. Ex- 
cept for the inner parts of the galaxies with long exposure 
times, these are strongly dominated by photon statistics and 
read-out noise. The contribution of additional sources of ran- 
dom error is quite negligible. 

(i) photon statistics and read-out noise: They have been 
computed with the help of a new set of analytic formulae 
derived by Cardiel et al. (1998). Error images, which are 
created at an early stage of the reduction procedure, are 
processed in parallel with data images. The basic arithmetic 
manipulations performed on the data frames are translated 
into the error frames following the usual error propagation 
laws. The errors obtained using this method are in excel- 
lent agreement with numerical simulations (see Cardiel et 
al. 1998 for details). In the outer parts of the galaxies the 
spectra were spatially binned to guarantee maximum ran- 
dom errors AD4000 ^ 0.15 and AMg2 Ss 0.025 mag. 

(ii) flux calibration: We observed four spectrophoto- 
metry standards (Massey et al. 1988; Oke 1990) in each 
run. The derived calibration curves of each run were aver- 
aged and we estimated the uncertainty in flux calibration as 
the rms scatter among the indices measured with different 
curves (AD4000 — 0.03, AMg2 ~ 0.004 mag for run 1, and 
AD4000 ^ 0.03, AMg 2 < 0.001 mag for run 2). 

(iii) wavelength calibration: Spectra were transformed 
to a linear wavelength scale by using ~ 35 arc lines fitted 
by 5th order polynomials, with typical rms ~ 0.3 A (run 1) 
and rms ~ 0.5 A (run 2). These scatters are completely con- 
sistent with those obtained by cross-correlating the central 
spectra of the single exposures of each galaxy. These un- 
certainties translate into AD4000 < 0.01 and AMg2 < 0.001 
mag for both runs. 

(iv) radial velocities: Radial velocities were computed 
by cross-correlating the central spectra of the final co-added 
image of each galaxy with different templates (usually ~ 4 
broadened spectra of bright G-K giant stars) . The different 
measurements were obtained with a typical rms ~ 10 km/sec 
(run 1) and rms ~ 20 km/sec (run 2), which has a negligible 
effect in the measurement of the indices (AD4000 < 0.01 and 
AMg2 < 0.001 mag). For comparison, an error of 200 km/sec 
turns into AD4000 — 0.01 and AMg2 ~ 0.004 mag for both 
runs. 



* Sec also: 

http: / / www.ucm.es / OTROS / Astrof/ rcduccme / rcduceme.htm] 



3.2 Systematic errors 

Systematic errors in the measurement of spectral indices 
arise mainly from the effects of spectral resolution (and 
velocity dispersion), flux calibration, sky subtraction and 
contamination by nebular emission, among others. Unfortu- 
nately, systematic errors do not allow the same straightfor- 
ward treatment used for the random errors. 

(i) spectral resolution and velocity dispersion: The Mg2 
index, which is defined using relatively broad bandpasses, 
is quite independent on spectral resolution and velocity 
dispersion corrections (e.g. Gorgas, Efstathiou & Aragon- 
Salamanca 1990; Gonzalez 1993; Carollo, Danziger & Buson 
1993). The D4000, due to its large wavelength coverage, is 
completely insensitive to these effects. Thus, these particu- 
lar indices do not require any spectral resolution correction. 

(ii) flux calibration: In order to guarantee that our Mg2 
measurements are in a consistent photometric system, we 
observed a sample of 39 (run 1) and 13 (run 2) template 
stars (typically G-K giant stars) from the Lick library (Gor- 
gas et al. 1993; Worthey et al. 1994). We compared our 
Mg2 indices with those of the Lick group in Fig. 1. We 
found a systematic offset of 0.021 mag (with a rms scatter of 
0.010 mag). This value is almost the same when analysing 
the data for both runs separately: 0.020 (rms = 0.008) and 
0.021 (rms = 0.021) mag for run 1 and 2, respectively. This 
offset is similar to that reported in CGA95 (0.014 mag), and 
it is also consistent with those determined by other authors: 
e.g. 0.017 mag (Davies, Sadler & Peletier 1993; Gonzalez 
1993), 0.010 mag (Carollo & Danziger 1994). This system- 
atic offset is due to the fact that the Lick data were flux 
calibrated using a tungsten lamp as a relative calibrator of 
the instrumental spectral response and, thus, they are not in 
a true photometric system. However, we have transformed 
our Mg2 indices to the Lick system to perform comparisons 
with Mg2 published in the literature. 

Unfortunately, the Lick library does not cover the spec- 
tral region corresponding to D4ooo- To check the reliability 
of the photometric system employed in the determination of 
this spectral feature we have compared the measurements 
of four stars in common with CGA95. The scarcity of this 
sample only allow us to conclude that the D4000 indices mea- 
sured here and those presented in CGA95 are, within errors, 
in the same photometric system (offset = 0.17, rms = 0.27). 

(iii) sky subtraction: Since galaxy light levels are usu- 
ally only a few per cent of the sky signal in the outer parts 
of the galaxies, this process constitutes one of the most im- 
portant potential sources of systematic errors when mea- 
suring line-strength gradients. This problem was treated in 
CGA95, and we refer the interested reader to that paper. 
Although there is not a simple recipe to detect a system- 
atic effect in the sky subtraction, anomalous deviations of 
the line-strength indices in the outer parts of the galaxies are 
indicative of its presence. We have corrected from the contri- 
bution of the galaxy light in the regions where the sky levels 
where measured by fitting de Vaucouleurs laws to the galaxy 
profiles. When necessary, we have subtracted from the sky 
spectra a scaled and averaged galaxy spectrum. Even in the 
outermost regions of the galaxies the effect of such correc- 
tions on the measured indices is nevertheless quite small 
(~ 10-20 per cent of the random error). 

(iv) contamination by nebular emission: Emission lines 
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present in the inner parts of some CCGs (e.g. [Nelll] A3869, 
[Sll] AA4069, 4076, H<5, [O m] A4959, [Ni] A5199) must 
be removed from the bandpasses employed in the mea- 
surement of Mg2 and D4ooo- For this purpose we have 
followed an interpolation procedure based on a simul- 
taneous fitting of a scaled template spectrum, a low 
order polynomial and a Gaussian, in the region close 
to the emission lines (see Fig. 2). The pixels affected 
by the presence of emission where carefully replaced 
by the result of that fit (i.e. polynomial+template). In 
some objects these corrections are not negligible: e.g. 
A 2597, <5D 4 ooo = 0.06, SMg 2 = 0.100; A 478, 5D 40 oo = 0.02, 
<5Mg 2 = 0.040; Hydra A, <5D 400 o = 0.01, 5Mg 2 = 0.010 
(<5 Index — Index nc i — Index e \, where Index c \ is the index 
measured without removing the emission lines, and Index nc \ 
the index obtained when the emission has been interpo- 
lated) . This procedure has been applied to the whole sample 
of galaxies with emission lines in their central regions. 



included in these diagrams. Filled circles and stars corre- 
spond to measurements in different sides of the galaxies. 
Open circles refer to secondary nuclei. The thick horizontal 
solid lines in the upper left corner of the panels indicate the 
spatial extent of the emission lines, when present. We define 
the emission line region as the area where the emission lines 
were detected above 1 a. The cut-off of this region is very 
sharp, so the actual deffinition of its edge is not critical. 

In most cases, indices measured at both sides of 
the galaxies agree within the errors. Error-weighted least- 
squares fits to all the points, excluding data within the cen- 
tral 1.5 arcsec (affected by seeing) and secondary nuclei, are 
shown as thin lines for all the galaxies. In addition, we have 
computed inner and outer gradients (i.e. inside and outside 
the emission line region, respectively), in the galaxies with 
central emission lines. These inner and outer gradients, plot- 
ted with a thick line, have been forced to join at the radius 
where emission lines end. The results obtained in all these 
fits, together with their formal errors, are given in Table 4. 



4 RESULTS 

4.1 Measurement of the central indices 

Following the same procedure than in CGA95, and in order 
to avoid aperture effects, we have obtained central D4000 and 
Mg2 indices employing a standard metric aperture size cor- 
responding to 4 arcsec projected at the distance of the Coma 
cluster (~ 2.6 kpc). The central indices presented in CGA95 
were not corrected from emission lines. In order to get a 
fully consistent list of central indices, we have revised the 
spectra of that sample which exhibited emission lines. We 
have found that the effect is non negligible in two objects, 
namely A 1795 (with 5D 4 ooo = 0.02 and <5Mg 2 = 0.026 mag, 
and A 2199 (with <5D 40 oo = 0.01 and <5Mg 2 = 0.009 mag). In 
addition, we have transformed all the Mg 2 measurements to 
the Lick system by applying the derived systematic offsets. 
A full list with all the available central measurements and 
associated random errors, including the revised data from 
CGA95, is given in Table 3. 



4.2 Line-strength gradients 

We have measured D 4 ooo and Mg 2 gradients for all the galax- 
ies of the sample. In order to obtain symmetric brightness 
profiles at both sides of the galaxy centre, the spectra were 
shifted by a fraction of pixel in the spatial direction. This 
process was not carried out with the central galaxy of A 2626 
since in this object a secondary nucleus was also observed 
inside the slit at only ~ 3.4 arcsec from the centre of the 
CCG. In this case a simultaneous fit of two Cauchy func- 
tions plus a low order polynomial allows us to determine 
the location of the centre of the galaxy and the secondary 
nucleus, together with the relative light contribution of each 
component nearby the central region of the galaxy. The spec- 
tra were spatially binned in the outer parts of the galaxies 
to guarantee a minimum value of the signal-to-noise ratio. 

D4000 and Mg 2 gradients are plotted in Figs. 3 and 4, 
respectively, where the galaxies are displayed in decreasing 
mass deposition rate (label in the upper right corner of each 
panel). Galaxies from CGA95 with reliable gradients are also 



4.3 Comparison with previous work 

We have carried out a comparison of our central D 4 ooo 
measurements in CCGs (see Table 3) with the available 
data in the literature. In Fig. 5 we plot D 4 ooo data from 
JFN87 and M089, and (U-b) nuc from M092 against our 
central D 4 ooo indices. There is a clear systematic offset of 
0.24 (rms = 0.04) in D 40 oo between the indices of JFN87 
and this work. Although error bars are not available for the 
JFN87 data, assuming a typical random error of AD4000 ~ 
0.04, the derived offset is probably inconsistent with the in- 
trinsic errors and must be considered as a systematic effect. 
On the contrary, the D 4 ooo values from M089 are consis- 
tent with our measurements since the computed offset is 
0.06 ± 0.08. The nuclear (U-b) colours from M092 dis- 
play a very good correlation with the D 4 ooo- A similar cor- 
relation was previously presented by M089. This correla- 
tion can be easily understood since the (U— b) colours mea- 
sure essentially the same spectral break than D 4 ooo- The 
dashed line is panel 5b is the least-squares fit to the data: 
(U-b)nuc = -1.03(±0.14) +0.81(±0.07)D4ooo- This linear 
relation will be employed later in this paper to estimate 
central D 4 ooo from (U— b) nuc colours. 

Line strength gradients of Mg in brightest cluster galax- 
ies have been presented by Fisher, Franx & Illingworth 
(1995, hereafter FFI95). Our full sample includes four galax- 
ies in common with theirs, namely N 2832 (A 779), N 6166 
(A 2199), N 7720 (A 2634, —note: this object was misla- 
belled as N 7728 in CGA95) and the central cluster galaxy 
of A 496. In order to compare their measurements of Mgb 
with ours of Mg 2 we have transformed Mgb into Mg 2 . The 
atomic Mgb index is very well correlated with the molec- 
ular Mg 2 index in early-type systems: Gonzalez & Gorgas 
(1995) report Mg 2 (mag) ~ 0.066 Mgb(EW). The compari- 
son with FFI95 data is shown in Fig. 6. The measurements 
for the central galaxy in A 779 exhibit an excellent agree- 
ment. However, this is not the case for the central galaxies 
of A 496, A 2199 and A 2634. 

The discrepancies in A 496 and A 2199 are readily un- 
derstood in terms of contamination by the emission line 
[N 1] A5199. Unfortunately, this line is centred in the red con- 
tinuum bandpass of the Mgb index, and its presence trans- 
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lates into an overestimation of the index value due to an 
enhancement of the measured continuum level. The effects 
of this emission have been studied in detail by Goudfrooij & 
Emsellem (1996). Since in Fig. 6 we have derived Mg2 from 
the Mgb indices of FFI95, these Mg2 values are also overes- 
timated. Although we cannot correct the FFI95 indices from 
the effect of the emission line, we can reproduce their mea- 
surements by measuring Mgb in our spectra without remov- 
ing the contamination of [Nl] A5199, and then transforming 
these Mgb values into Mg2 indices. The effect of such pro- 
cedure is shown in panels 6b and 6c as dotted arrows. It is 
clear from the figures that in this case the agreement be- 
tween FFI95 data and our gradients is then fairly good for 
A 2199 and A 496. In fact our indices lie between those 
determined by FFI95 along the major and minor axis. 

It is important to note that the effect of [N i] A5199 in 
Mg2 is much smaller than in Mgb. Goudfrooij & Emsellem 

(1996) showed that a strong EW of [N i] 5199 ~ 1.8 A trans- 
lates into AMgb ~ 2 but only AMg2 ~ 0.03 (transforming 
this AMgb into AMg2 using the linear relation described 
above gives AMg2 ~ 0.13, which is four times larger). 

Nevertheless, the Mg2 gradient for the central galaxy 
of A 2634 does not show the same agreement found in the 
previous objects (note that, in this case, Mg2 is not affected 
by the emission lines). Our measurements are clearly below 
the FFI95 data. Additional central Mg2 measurements for 
A 779, A 2199 and A 2634 are also available from Trager 

(1997) — plotted as a short horizontal full line — and Lucey 
et al. (1997) — short horizontal dashed line. These central 
indices agree with our central values for A 779 and A 2199, 
confirming the previous comparison with the FFI95 data. 
In addition, the central Mg2 values of A 2634 from Trager 
(1997) and Lucey et al. (1997) he in the region between the 
FFI95 data and ours. This discrepancy, which affects mainly 
the absolute values of the Mg2 index, but not the shape 
of the radial gradient, is not understood, so the caution is 
needed when interpreting the A 2634 data. 



5 DISCUSSION 

5.1 Central spectral indices 

Different authors have reported correlations of central blue 
colours and spectral features with mass deposition rates: 
JFN87 found a correlation between D4000 and the mass flow 
rate (M v ) within their spectrograph slit; M089, M092 and 
McNamara (1997) shown that central U — B colour excesses 
in CCGs immersed in clusters with high mass deposition 
rates also correlate with M; in CGA95 we confirmed the 
correlation with D4000 and presented a correlation with Mg2. 

In Fig. 7 we have plotted all the central measurements 
listed in Table 3, as a function of the X-ray derived mass 
deposition rate. In panels 7a and 7b we have used the to- 
tal M values, whereas in panels 7c and 7d we computed 
the central mass deposition rates in the same fixed metric 
aperture size employed in the measurement of the central in- 
dices, assuming M oc r. We have also included in panels 7a 
and 7c central D4000 estimates of eight additional galaxies 
obtained from the nuclear (U— b) colours of M092. For this 
purpose we have used the linear fit derived in Section 4.3 



These new indices are listed in Table 5. Filled and open sym- 



bols in Fig. 7 correspond to CCGs with and without central 
emission lines, respectively. 

Although, as mentioned above, central D4000 and Mg2 
indices are well correlated with mass deposition rate in the 
sense that both indices decrease when the mass flow in- 
creases, it is interesting to note that some objects lie off the 
main relation. In particular the cooling flow galaxies without 
emission lines A 644, A 2029 and A 2142 are found above 
the relation, being their central indices fully consistent with 
those of central galaxies in clusters without cooling flow. Si- 
multaneously, there are also two cooling flow galaxies with 
emission lines, Hydra A and N 1275, which also seem to 
depart from the relation with M. 

To perform a simultaneous study of the variations in 
D4000 and Mg2, we present an index-index diagram in Fig. 8. 
In panel 8a we plot the central indices with their errors for 
the whole sample. In the following panels of this figure we 
have suppressed the galaxy names and error bars for clarity. 
Since the CCGs exhibit a clear sequence in this index-index 
plane, in panel 8b we compare this sequence with that dis- 
played by early-type galaxies. The small dots correspond to 
a sample of 135 elliptical and SO galaxies observed by Kim- 
ble, Davidsen and Sandage (1989; note: we have transformed 
their Mgb dat a into Mg2 using the same relation employed 
in section 4.3). As we already shown in CGA95, CCGs in 
cooling flow clusters populate the D4ooo _ Mg2 plane following 
a different sequence to that displayed by elliptical galaxies. 
The enlargement of the galaxy sample allow us to gain new 
insights into this study. In particular, it is clear from the 
figure that the CCGs in cluster without cooling flows (open 
circles) lie in the same region than the giant ellipticals (high 
Mg2 and D4000 values) . The central galaxies in clusters with 
cooling flows which do not exhibit emission lines in their 
central parts (open stars, with symbol size proportional to 
mass deposition rate) are also compatible with the location 
of the non cooling flow galaxies, even for high M, like A 2142 
and A 2029. On the other hand, cooling flow galaxies with 
central emission lines (filled circles, with symbol size pro- 
portional to M) depart from this location, being the largest 
accretors the galaxies which show the largest deviations. Us- 
ing the same symbols, we plot in panel 8c the D4000 obtained 
from (U— b)nuc colours of M092 (as explained above). It is 
interesting to note that some cooling flow galaxies with only 
moderately high M (Hydra A and N 1275) display quite 
low D4000 values. Thus, it is clear that cooling flow galaxies 
with emission lines do not represent a simple one-parameter 
sequence in mass deposition rate. 

We investigate now the behaviour of the spectral indices 
as a function of age and metallicity by studying the evolu- 
tion of single bursts of star formation. We have used the 
galaxy isochrone synthesis spectral evolution library, GIS- 
SEL96, from Bruzual & Chariot (1997). In panel 8d we plot 
the evolution of three single bursts of star formation (full 
lines) with different metallicities ([Fc/H] =— 0.64, —0.33 and 
0.09). The model sequences are labelled with small circles 
which indicate the age of the burst in Gyr. Mg2 values in 
the models are computed using the empirical fitting func- 
tions derived by Gorgas et al. (1993) and Worthey et al. 
(1994). Unfortunately, this type of fitting functions are not 
available for the D4000 , and a non- negligible systematic offset 
between the predictions of the models and the data could be 
present. For simplicity, we represent the sequence of normal 
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ellipticals and SO of Kimble et al. (1989) by a dashed par- 
allelogram. Taking into account the models, we have drawn 
two vectors to indicate changes of age and metallicity in this 
D4000 - Mg2 plane and, interestingly, both parameters are not 
completely degenerated in this diagram. It suggests that nor- 
mal ellipticals represent a metallicity sequence, whereas the 
cooling flow objects describe a locus quite parallel to the age 
evolution. This strongly suggests that the observed spectral 
properties of the cooling flow galaxies can be understood as 
a consequence of star formation. 

Following the procedure described in Gorgas et al. 
(1990) and in CGA95, we have modelled the expected 
changes in the central D4000 and Mg2 indices when a nor- 
mal CCG, initially without cooling flow, undergoes star for- 
mation as a result of the cooling flow — panel 8e. All the 
models describe the effect of such star formation overim- 
posed on a galaxy with D4000 and Mg2 indices equal to the 
averaged values obtained with our subsample of five CCGs 
in clusters without cooling flows (i.e., typical of old stellar 
populations). Dotted lines represent the effect of a contin- 
uous star formation during the last 0.1 Gyr, 1.0 Gyr and 
10.0 Gyr (from bottom to top). Full lines correspond to sin- 
gle bursts occurred 0.01 Gyr, 0.1 Gyr and 1.0 Gyr ago (also 
from bottom to top). In all these models a Scalo (1986) IMF 
and solar metallicity were adopted. Finally, the dashed line 
indicates the effect of a power- law spectrum (Fx <x A Q_2 , 
with a = 0.5). No intrinsic reddening has been taken into 
account in the models. Small dots in the model lines indi- 
cate the fixed fraction fy (given by the numbers) of the 
V light that comes from the accretion population (i.e., the 
newly-formed stars) . In the case of a continuous star forma- 
tion, this fraction fy can be related to the SFR using the 
expression 



SFR = 



fv Ly (M/L)f 
t 



(1) 



where Ly is the total V luminosity of the galaxy, (M/L) ap 
is the stellar mass-to-light ratio of the accretion population, 
given by the models, and t is the time over which the star 
formation has taken place. Similarly, the parametrization for 
the models involving a single burst of star formation is given 
by 



Map = fv L v (M/L), 



(2) 



where Map is the total mass in stars formed in the accretion 
population. 

Panel 8e reveals that the sequence of cooling flow galax- 
ies does not seem to be very well reproduced with a con- 
tinuous star formation during the last 1-10 Gyr. Either a 
young (< 1.0 Gyr) continuous star formation, or simply a 
recent burst of star formation ~ 0.1 Gyr ago, can match 
the observed central indices. In the lower right corner of 
panel 8e we have also plotted the effect of intrinsic redden- 
ing in the measurements. As an example, we represent the 
corrections for an hypothetical object with D4000 = 1-7 and 
Mg 2 = 0.300 as a function of E(B - V) (the colour ex- 
cess ranges from 0.0 to 1.0, with ticks spaced 0.1 units). As 
it could be expected, the effect on Mg2 is almost negligible, 
whereas the contrary is true for the D4000 ■ Although Hansen 
et al. (1995) have reported E(B -V)< 0.1 for Hydra A, 
other authors have found non-negligible reddening values 
(E(B — V) ~ 0.1-0.6) for several cooling flow galaxies (e.g. 



Hu 1992; Allen et al. 1995; Allen 1995; Donahue & Voit 1993; 
Voit & Donahue 1997). If we allow for the existence of inter- 
nal reddening, the case for recent bursts instead of a continu- 
ous star formation is even stronger. Nevertheless, it is impor- 
tant to keep in mind that the effect of reddening in the stellar 
continuum is usually lower than in the emission lines, e.g. 
Calzetti (1997) reports E(B - V) s tellar = 0.44 E(B - V) gas 
in typical star forming regions. 

However, we need to be cautious since, apart from in- 
trinsic reddening, there are additional sources of uncertainty 
present in the models. First, there is some uncertainty in the 
starting point of the model sequences, since the dispersion in 
the central indices of CCGs without cooling flows is probably 
real. Second, we have already mentioned that the modelling 
of D4000 is somewhat uncertain given the lack of empirical 
calibrations, thus some systematic effects could be present. 
And finally, the models assume solar metallicity and solar 
element ratios: if there is an enhancement of [Mg/Fe] in the 
cooling flow material, the predicted Mg2 indices could be un- 
derestimated. For these reasons, it seems premature to rule 
out competely, from the models shown in panel 8e alone, a 
continuous star formation (lasting some Gyrs) as the source 
of the blue light. 

Considering the relatively short wavelength baseline be- 
tween the D4000 and Mg2 indices, the contribution of a 
power-law spectrum could also reproduce the observed spec- 
tral changes. This result is not surprising since, in a short 
wavelength interval, the spectrum of an early-type star can 
be approximated quite well by a power law (see also Craw- 
ford & Fabian 1993). However, as we have already discussed 
in the introduction, the absence of polarization in the blue 
light of A 1795 (McNamara et al. 1996b), does not give 
strong support to this alternative. In addition, some cen- 
tral galaxies show clear signatures of strong Balmer ab- 
sorption lines (e.g. N 1275, Crawford & Fabian 1993; Hy- 
dra A, Hansen et al. 1995, Melnick et al. 1997). Allen (1995) 
found that the strong UV/blue continua in the spectra of 
a small sample of CCGs were better described by young 
stars than by power-law emission models. This author also 
reported positive detection of Wolf-Rayet features in two 
CCGs (those in A 1068 and A 1835), which favour the idea 
of the formation of massive stars. In this sense, Voit & Don- 
ahue (1997) have measured a weak Hell A4686 emission in 
A 2597. We also find evidences of such helium feature in the 
central spectrum of the dominant galaxies in A 1795 and in 
A 2597 (see Fig. 9) 

The role of metallicity and the choice of IMF on the 
model predictions plotted in panel 8e is explored in panel 8f. 
We show the GISSEL96 (Bruzual & Chariot 1997) predic- 
tions for different metallicities (as explained in the panel leg- 
end), and two IMFs, namely Salpeter (1955) — thin lines — 
and Scalo (1986) — thick lines — , for the single burst model 
with an age of 0.1 Gyr. The predictions of the models do not 
depend strongly on the choice of both [Fe/H] or the IMF. 

In order to obtain further constraints on the range 
of acceptable models, we need to explore a larger wave- 
length range. We have done so using UV and near-IR data 
from the literature (see Fig. 10 and description in the fig- 
ure caption). Although a continuous star formation last- 
ing 10 Gyr could explain the observed optical and near- 
IR observations, the UV fluxes are underpredicted. Recent 
bursts with ages ~ 0.1-0.2 Gyr are simultaneously compat- 
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ible with both the UV and IR fluxes, and with the mea- 
sured central line-strength indices. Note however, that at 
UV wavelengths the effects of dust extinction can be quite 
dramatic even for moderate colour excesses: using the in- 
terstellar reddening curve of Savage & Mathis (1979), and 
assuming E(B — V) = 0.2 and E(B - V) = 0.4 we estimate 
that the measured UV fluxes need to be corrected upwards 
by factors of ~ 5-6 and ~ 30-20, respectively. In these cases, 
is quite obvious that the UV continuum would constrain the 
upper limit of the burst age to be < 0.1 Gyr. This result 
agrees with the behaviour of D4000 and Mg2 with E(B — V) 
in Fig 8e, indicating that a recent burst would be strongly 
favoured against constant star formation. 

Another piece of evidence which can help to understand 
the behaviour of the spectral indices of the CCGs comes 
from the emission line properties. Using line-ratio diagrams 
HBvBM classified the emission nebulae at the centres of cen- 
tral dominant galaxies into two types. In particular, Type I 
corresponds to objects with relatively high [N 11] A6584/Ha 
(~ 2.0) and [S Il]A6717/Ha (~ 0.7), whereas Type II ob- 
jects exhibit lower line-ratios ([N ii]A6584/Hq ~ 0.9 and 
[S n]A6717/Ha ~ 0.4). In panels 8g and 8h we have plotted 
the available HBvBM types for the galaxies in our sample 
and in M092 respectively. There is a clear dichotomy in the 
way Type I (plotted with "1" symbols) and Type II (plot- 
ted with "2" symbols) galaxies populate the central cluster 
galaxy sequence. This result is not surprising since Type I 
objects are associated to clusters with low X-ray and Ha 
luminosities, small optical nebulae and low mass deposition 
rate, whereas the contrary is true for Type II objects (HB- 
vBM, Baum 1992, Donahue & Voit 1997). 

Another clue comes from the morphology of the blue 
excess light. Recently, McNamara (1997) has established a 
morphological classification based on the study of the sur- 
face brightness distributions and colour structure in cooling 
flow CCGs. This classification consists in four morphologi- 
cal types ordered by increasing geometrical complexity: (1) 
point source (e.g. A 2199, A 2052) with unresolved blue 
nuclei and nebular line emission, low nebular luminosity, 
with modest blue anomalies; (2) disk (e.g. Hydra A) un- 
usual in cooling flow galaxies and characterised by a disk of 
gas and young stars rotating around the nucleus; (3) lobe 
(e.g. A 1795, A 2597) with bright, blue lobes of optical con- 
tinuum separated several kiloparsecs from the nucleus; and 
(4) irregular-amorphous (e.g. N 1275 may be in transition 
between types 3 and 4). In Fig. 8g we have also plotted the 
four morphological types in the neighbourhood of the galax- 
ies which display such characteristics. Taking into account 
the D4000 exhibited by N 1275, we have tentatively located 
Type 4 near PKS0745 - 191. Very interestingly, McNamara 
(1997) has speculated about the possibility that the struc- 
tural types reflect actually an evolutionary sequence. We will 
discuss this point in more detail in section ^| 

5.2 Line-strength gradients 

Although line-strength gradients in Mg2 (e.g. Faber 1977; 
Gorgas et al. 1990; Gonzalez 1993; Davies et al. 1993; Carollo 
et al. 1993; FFI95; Gonzalez & Gorgas 1995) and in D4000 
(e.g. Munn 1992; Davidge & Grinder 1995) are usually found 
to be linear with log(r) in elliptical galaxies, in CGA95 we 
showed that some cooling flow galaxies exhibit a clear slope 



change at intermediate radii. With the inclusion of the new 
data presented in this paper (see Figs. 3 and 4), we find 
that this change of slope is clearly related to the presence of 
emission lines in the central regions of CCGs. 

We measured the radius at which the change of slope 
takes place using error-weighted least-squares fits to two 
straight lines forced to join at a variable radius. The method 
gives the break radius which minimises the residual variance 
of the fits. We estimated the uncertainty in the break radius 
using Monte Carlo simulations that take into account the 
measurement errors in the indices (typically ~ 1000 sim- 
ulations per gradient). For the emission line galaxies with 
reasonably good data beyond the emission region, we have 
compared in Fig. 11 the break radius with the observed 
emission-line radius (as defined in section 4.2). The figure 
reveals that the break radii are almost coincident with the 
emission-line radii, with the former being slightly larger. 

It is quite evident from Figs. 3 and 4 (see also the fitted 
gradients tabulated in Table 4), that cooling flow galaxies 
with emission lines exhibit flat and even positive gradients 
in the region where the emission is detected (e.g. A 478, 
A 1795, A 2597, Hydra A, A 85, A 2199). In fact, these in- 
ner gradients (measured in the emission line region) seem 
to be correlated with the mass deposition rate, i.e. more 
positive with increasing M (see Figs. 12a and 12b). In the 
outer parts of these galaxies, where emission is not observed, 
the derived mean D4000 gradient is consistent with that ob- 
served in normal elliptical galaxies (see Fig. 12c). However, 
the mean Mg2 gradient in the outer regions of cooling flow 
galaxies with emission lines seems to be marginally stepper 
than the mean gradient for normal ellipticals (Fig. 12d), al- 
though the quality of the Mg2 data at large galactocentric 
distances is not good enough to reach a firm conclusion. A 
couple of galaxies (A 1126 and A 262) seem to have slightly 
stepper gradients in the inner parts than in the outer re- 
gions, although this result is not very significant due to the 
poor S/N in the outer spectra. In addition, mean gradients 
in central dominant galaxies with cooling flow but without 
emission lines (e.g. A 644, A 2142) are similar to those in 
galaxies without cooling flows and in elliptical galaxies. 

If we assume that the observed slope changes in the line- 
strength gradients are due to radial index variations induced 
by star formation in the emission-line region, we can use 
that information, coupled with the star formation models, 
to estimate the amount and spatial profile of the accretion 
population. We assume that the gradients determined in the 
outer regions represent the underlying gradient, i.e., that of 
the galaxy before the cooling-flow induced star formation 
takes place. That is justified by the similarity of these outer 
gradients with those of "quiescent" galaxies (i.e., both nor- 
mal giant ellipticals and CCGs without cooling flows). In 
this analysis we have exclusively concentrated on the D4000 
gradients, which are much better determined than the Mg2 
ones. 

As discussed in section 5.1, a single burst model with 
an age ~ 0.1 Gyr, solar metallicity and Scalo IMF is able 
to reproduce the observed spectral properties of the cool- 
ing flow galaxy sequence. We have employed this particular 
model to tentatively obtain the spatial profile of the accre- 
tion population. First we computed the averaged indices in- 
side the emission line region (simulating a circular aperture 
by weighting with light and radius), both from the fitted 
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inner gradient and from the extrapolation of the outer gra- 
dient into the emission line region. Then, we used the star 
formation model to derive the mean fv of the accretion pop- 
ulation required to reproduce the observed changes in the 
indices (i.e., the difference between the measured values and 
those estimated from the extrapolation of the outer gradi- 
ent). Assuming that the absolute magnitude corresponding 
to the old stellar population in the emission line region is 
(My) ~ —22.3 magj^, the mean fy value is then employed 
to compute the total V luminosity (Ly) of the galaxy plus 
the burst. Using Eq. |^, and the mass-to-light ratio of the ac- 
cretion population given by the model, we obtain the total 
mass in stars of the accretion population formed in the last 
burst. Note that, since the excess blue light is completely 
dominated by the youngest stars, this method is only sensi- 
tive to the stars formed in the most recent burst. 

In Fig. 13 we have plotted the estimated total amount 
of mass transformed into stars in the emission-line region 
versus the mass deposition rate computed in the same area 
(assuming that M ocr). There is an increase in the mass 
transformed into stars with M, with Hydra A lying some- 
what outside this trend. Interestingly, the total mass ac- 
creted by the cooling flow in a ~ 0.1 Gyr period inside the 
emission line region is remarkably similar to the mass in new 
stars (see section 

Assuming spherical geometry, it is possible to apply the 
same technique to obtain the amount of mass transformed 
into stars in concentric shells at different radii. Using a de- 
projection algorithm described in Appendix A, we have de- 
rived the spatial profiles for the density of mass transformed 
into stars (Fig. 14). We also show (dashed line) the expected 
density profile of the mass deposition assuming M(< r) oc r. 
Interestingly, the derived profiles seem to be quite parallel to 
the dashed line (the mean slope for all the galaxies excluding 
Hydra A is —2.08, rms = 0.18), whereas Hydra A exhibits 
a statistically significant steeper trend. This result strongly 
suggests a direct connection between the mass deposition 
and the star formation profiles. Indeed, figures 13 and 14 
suggest a very close link between the star formation and the 
cooling flow phenomenon. Such close connection would be 
difficult to understand without most of the gas accreted in 
the emission-line region being transformed into stars. 



t Hoessel, Gunn & Thuan (1980) reported that the mean mag- 
nitude of a sample of brightest cluster galaxies inside a metric 
radius of 16//ieo kpc is {My) = —22.68 ± 0.03 mag. In addi- 
tion, Hoessel (1980) showed that the internal regions of brightest 
cluster galaxies are adequately described by a modified Hubble 
law of the form I(r) = / c /(l + r' 2 /r 2 ), where I c is the central 
intensity, r c the core radius and r the angular distance from 
the centre. Assuming this intensity profile, the integrated lumi- 
nosity inside a radius r is then given by L(r) = it I c r 2 ln(l + 
r 2 /r 2 ). Introducing the value of the Hubble constant adopted 
in this paper (H = 50 km s — 1 Mpc - 1 ), and the mean core ra- 
dius derived by Hoessel (1980), (r c ) = 2.23//i60 kpc, we obtain 
(My) ~ —22.3 mag inside a radius r ~ XX/hso kpc, which is the 
average radius of the emission line regions in our sample. 



6 SUMMARY: TOWARDS A COMPLETE 
PICTURE? 

Many pieces of the cooling flow puzzle are now available. 
We summarise some of them here before trying to sort out 
a scenario able to explain the observational data. 

(i) Central D4000 and Mg2 indices in cooling flow galax- 
ies are correlated with mass deposition rate when emission 
lines are found in the central regions of such galaxies. 

(ii) The nuclear indices of the cooling flow galaxies in 
our sample which do not exhibit emission lines do not fol- 
low this correlation with M. The central line-strengths of 
both CCGs in clusters without cooling flows and cooling 
flow galaxies without emission lines are consistent with the 
values observed in giant ellipticals. 

(iii) The central D4000 and Mg2 indices of elliptical 
galaxies define a relatively narrow trend in an index-index 
plane. Cooling flow galaxies with emission lines exhibit a 
clear sequence in this diagram which departs from the locus 
defined by elliptical galaxies. When interpreted using stel- 
lar population models, this provides strong evidence for star 
formation in these galaxies. 

(iv) Although line-strength gradients are usually found 
to be linear with log(r) in elliptical galaxies, D4000 and Mg2 
gradients in cooling flow galaxies with emission lines show 
a clear slope change in the region where the emission is de- 
tected. In fact, these inner gradients seem to be correlated 
with the mass deposition rate. In the outer parts of these 
objects, where emission is not observed, the derived mean 
line-strength gradients are consistent with those in elliptical 
galaxies. 

(v) Mean gradients in cooling flow galaxies without 
emission lines are similar to those in CCGs without cool- 
ing flow and those in giant elliptical galaxies. 

(vi) The presence of emission lines must be in some way 
related to the cooling flow phenomenon since H/3 (JFN87) 
and Ha luminosities (HBvBM; Donahue & Voit 1997) are 
correlated with mass deposition rate. However, considering 
the high scatter in line luminosities for a given M, the rela- 
tionship does not seem to be a simple one (Baum 1992). 

(vii) Whichever mechanism is responsible for the di- 
lution of the measured spectral indices, the correlations of 
the central indices and gradients with the mass deposition 
rate in the emission-line region add further weight to the 
idea of a link between the blue excess and the cooling flow. 
Nevertheless, as Crawford & Fabian (1993) pointed out, the 
fact that there are central galaxies in clusters with high M 
but without blue excess suggests that the cooling flow alone 
cannot be responsible for the blue light. In particular, these 
authors noted the existence of strong cooling flow galaxies, 
with (e.g. A 478) and without (e.g. A 2029) emission-line 
nebulae, which apparently did not exhibit a blue continuum 
excess. However, our D4000 and Mg2 gradients in the cen- 
tral galaxy of A 478 does exhibit a clear positive behaviour 
(i.e. blue excess) in the emission-line region. An important 
question is whether there are any cooling flow galaxies with 
emission-line nebulae which do not exhibit a blue excess, 
and vice versa. So far this does not seem to be the case. 

(viii) Through a detailed analysis of the emission-line 
nebula in the central cluster galaxy of A 2597, Voit & Don- 
ahue (1997) conclude that hot stars are the most likely ion- 
izing source of the gas if some mechanical form of heating 
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(flowing from the hot intracluster medium into the cooler 
nebula) supplements photoelectric heating. 

(ix) With the available data, there is a clear dichotomy 
in the way cooling flow galaxies populate the D4ooo~Mg2 
plane depending on their emission-line properties as de- 
scribed by the classification scheme of HBvBM. The relative 
position of Type 1 and Type If objects in Fig. 8g strongly 
suggests that recent star formation (and therefore photoion- 
ization by hot stars) is the key to understand the observed 
dichotomy. 

(x) The excess of blue light (i.e. line-strength dilution) 
in the continuum spectra of CCGs is consistent with recent 
star formation episodes (e.g. Bertola et al. 1986; JFN87; 
Shields & Filippenko 1990; Crawford & Fabian 1993; Hansen 
et al. 1995; Allen 1995; Melnick et al. 1997; Figs. 8e and 10 
— this work). 

(xi) If the observed index variations are due to star for- 
mation in the emission-line region, the existence of correla- 
tion between the central indices and the line-strength gradi- 
ents with the mass deposition rates implies that the cool gas 
must be transformed into stars with a similar efficiency from 
galaxy to galaxy. In addition, the radial profiles of the star 
formation density (Fig. 14) are remarkably parallel to the 
expected density of the mass deposition profile. Both results 
are difficult to understand unless a large fraction (probably 
most) of the cooling flow gas accreted into the emission-line 
region is converted into stars. 

(xii) A large fraction (60-70%) of cD galaxies in cool- 
ing flows are radio- loud, whereas this probability decreases 
to 20% for cD galaxies in non-cooling flow clusters, and 14% 
for typical elliptical galaxies in clusters (Burns et al. 1997). 
Since the kinetic energy of the radio plasma can be a sig- 
nificant fraction of the thermal energy in the cooling flow, 
the correlation between radio emission and the presence of 
cooling flows is probably indicating the existence of inter- 
actions between the radio plasma and the centres of such 
cooling flows. However it is important to note that there is 
only a weak correlation between the 1.4 GHz radio power 
(HBvBM) and the 6-cm radio power (Burns 1990) with M. 

(xiii) The CCGs in A 1795 and A 2597 host blue lobes 
of optical continuum along the edges of their radio lobes 
(McNamara & O'Connell 1993; Sarazin et al. 1995; McNa- 
mara et al. 1996a, 1996b) . Through the analysis of numerical 
simulations De Young (1995) has shown that the blue con- 
tinuum can be due to young stars produced by jet-induced 
star formation. This author has suggested that the most ap- 
pealing scenario consists in short-lived jets (with ages less 
than 10 7 yr), with the high-pressure cooling flow being re- 
sponsible for the disruption of the jets. 

(xiv) Using an X-ray colour/deprojection technique, 
Allen & Fabian (1997a, 1997b) have proved that all the cool- 
ing flow galaxies in their sample exhibit significant central 
concentrations of cooling gas. These authors find that in- 
trinsic X-ray absorption increases with decreasing radius, 
and that significant excess absorption is not present in non 
cooling-flow clusters. The estimated times required for the 
cooling flows to accumulate the observed absorbing gas are 
typically of only a few 10 s years. 

Summarizing, all these results suggest that very likely 
the presence of cooling flows, emission lines, blue excesses, 
radio emission, soft X-ray absorbing gas and star forma- 
tion are different aspects of the same phenomenon. As we 



have previously mentioned (section p.l[ Fig. 8g), McNamara 
(1997) has suggested that his morphological classification of 
the blue structure in cooling flow CCGs reflects an evolu- 
tionary sequence. This author proposes the following sce- 
nario: central cluster galaxies accrete ~ 10 s Mq of gas from 
a gas-rich dwarf galaxy or a cooling parcel of hot gas. The 
accreted gas flows to the galaxy centre where it triggers a 
radio source (by providing the necessary fuel), which subse- 
quently induces a burst of star formation due to the rapid 
collapse of cold clouds compressed by shocks along the ex- 
panding radio source (type 3, lobe). In a few ~ 10 7 yr the 
lobes disperse (type 4, amorphous), leading finally to an un- 
resolved point source (type 1). The type 2 (disk) objects 
do not fall inside this evolutionary sequence. As McNamara 
(1997) points out, Hydra A is an example of tidal accre- 
tion and likely not accretion from the cooling flow. In this 
sense, Melnick et al. (1997) have proven that the nuclear 
emission lines in this galaxy originate in a gaseous disk-like 
structure with a rotational velocity of ~ 300kms _1 . Under 
this hypothesis, it is not surprising that the behaviour of 
the central galaxy of Hydra A in Figs. 7, 12, 13 and 14 is 
different from that of the remaining cooling flow objects. 
The locus of N 1275 in Figs. 7a and 7c is also suspicious. 
However this object exhibits such unusual properties that it 
has been interpreted either as a system of colliding galaxies 
(Baade & Minkowski 1954; Minkowski 1957; Rubin et al. 
1977; Hu et al. 1983; Holtzman et al. 1992), two superposed 
galaxies without interaction (De Young, Roberts & Saslaw 
1973), or even an exploding galaxy (Burbidge, Burbidge & 
Sandage 1963; Burbidge & Burbidge 1965). Thus, it is likely 
that the morphology displayed by this galaxy is more com- 
plicated than that naively expected from the cooling flow 
phenomenon. 

If this scenario is true, an important question is why 
are we observing so many objects with signatures of recent 
star formation. Given the relatively large fraction of type 3 
objects observed in our sample (at least two: A 1795 and 
A 2597), and considering the timescale for the transition 
of a CCG throughout this stage (a few ~ 10 7 yr), the sce- 
nario proposed by McNamara (1997) is not tenable unless 
radio triggered star formation takes place several times dur- 
ing the lifetime of the central galaxies. Some indication that 
this could well be the case comes from the study carried out 
by Allen (1995) who found that the number ratios of stars 
of different spectral types needed to reproduce the UV/blue 
continua in a sample of CCGs also favour a scenario in which 
star formation occurs in a series of bursts. The analysis we 
have presented here (section 5) also suggests that a young 
burst (~0.1 Gyr of age) is able to explain the observed spec- 
tral properties of the emission-line CCGs better than a more 
continuous star formation. The population synthesis models, 
in conjunction with the spectral data, also indicate that the 
total mass of cooled gas transformed into stars is roughly 
equal to the mass accreted inside the emission-line-region 
over a ~ 10 8 yr period. Interestingly, the accumulation time- 
scales derived by Allen & Fabian (1997a, 1997b) to repro- 
duce their calculated masses of absorbing gas is also only a 
few 10 s yr. Since the observed blue light would be completely 
dominated by the stars formed in the last burst, it is possible 
that a sequence of bursts has occured before the most recent 
one, fueled by the accreted gas. Since cooling flows can be 
disrupted by cluster-cluster interactions (Burns et al. 1997), 
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it is also possible that the sequence of bursts takes place 
only occasionally. As Crawford (1997) pointed out, if rich 
clusters were built up hierarchically, the subsequent merg- 
ers of the subclusters could have disrupted earlier cooling 
flows, likely high-redshift radio sources (Bremer, Fabian & 
Crawford 1997). Therefore, strong cooling flows could have 
existed in today's rich clusters only over the last 10 9 years, 
allowing for a few bursts of star formation (spaced by a few 
10 s years) having taken place during this period. 

This scenario accounts for the fate of the gas accreted 
within the region where emission lines are observed (i.e., 
the inner ~ 20kpc): it accumulates over periods of a few 
10 8 years and then most of it (perhaps all) forms stars in 
bursts. Outside this region, star formation signatures are 
not observed, perhaps because the radio source is unable to 
trigger star formation at large distances. What happens to 
the gas depossited there is not clear: some of it is observed 
as cool material out to ~ 100 kpc (Allen & Fabian 1997a, 
1997b), but our lack of knowledge of the physics of star 
formation prevents us from understanding why (massive) 
stars are not formed. 

This scenario could also provide a tentative explanation 
of why there are CCGs in clusters with strong cooling flows 
that do not show any signs of recent star formation: a pow- 
erful radio source is also necessary to trigger or catalise the 
formation of (massive) stars. The cooling flow provides the 
fuel, but without a central engine the radio source, needed to 
kick-start the star formation, would not develop. The poor 
correlation between radio power and star formation indica- 
tors (e.g., Ha fluxes) could be accounted for by the difference 
in time-scales involved in the radio phenomenon (~ 10 7 yrs) 
and the star-formation signatures (~ 10 8 yrs). There could 
also be some delay between the onset of a powerful radio 
source and the begining of the star formation, which could 
explain the few cases where a radio source is present in cool- 
ing galaxies without clear signs of excess blue light. 

We must stress that this picture, although very sugges- 
tive, is highly speculative, and further observational work 
is clearly needed. Our study of optical absorption features 
has proven to be an excellent tool to investigate the effect of 
cooling flows in central dominant galaxies, but we recognise 
that extending the wavelengh baseline of the observations is 
necessary to provide stronger constraints. Moreover, a sys- 
tematic study of the spatial extent and morphology of the 
star formation (e.g., via Ha images) and its correlation with 
radio emission would provide essential information on the 
suggested link between radio activity and star-formation. 
Near-infrared images would add to the picture by charac- 
terising the distribution of the underlying old stars. 

Progress is also required in the detailed modelling of 
the behaviour of the spectral features for potentially com- 
plex stellar population with a mixture of ages, metal abun- 
dances and element ratios. A missing ingredient is a reliable 
empirical calibration of some spectral features with intrinsic 
stellar parameters. At present we are near to complete the 
D4000 calibration in the Lick stellar library (Cardiel, Gorgas 
& Pedraz 1998). 
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APPENDIX A: RADIAL PROFILE OF THE 
STAR FORMATION 

We derive the spatial profile of the density of the mass trans- 
formed into stars (Fig. 14) assuming that the change of slope 
observed in the D4000 gradients in CCGs with cooling flows 
and emissio n lin es are due, exclusively, to star formation 
(see Section 5.2). In order to avoid uncertainties in the ab- 



solute values of the total V luminosity of the CCGs in the 
emission-line regions, we only derive the relative radial varia- 
tions (normalizing both the density of the mass transformed 
into stars and the spatial extent of the emission-line region). 
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The extrapolation of the outer D4000 gradients into the 
emission-line region allows us to estimate the expected ra- 
dial behaviour of the D4000 corresponding to the underlying 
old stellar population. Assuming a particular star formation 
model (in this case we have chosen a single burst with an 
age ~ 0.1 Gyr, with solar metallicity and Scalo IMF), it is 
possible to derive the fv required to reproduce the observed 
D4000 decrease (i.e. the difference between the measured and 
the extrapolated D4000 values) as a function of radius. 

If we assume that galaxies exhibit spherical geometry, 
for any concentric shell Eq. can be written as 

M AP (k) = f v (k) L v (k) (M/L) A p, (Al) 

where k is an integer which indicates the number of the shell 
(numbered from the center to outside), M\p{k) and fv(k) 
are the total mass transformed into stars and the fraction 
of V light that comes from the accretion population in the 
fc th -shell (respectively), and (M/L)ap the mass-to-light ra- 
tio of the accretion population given by the adopted stellar 
population model. The total V luminosity of the fc th -shell (in 
arbitrary units) is computed from the number of counts in 
the V wavelength range within the spectrograph slit, Ny(k), 
at a radius r(fc) (in arbitrary units) using 

L v {k) = Nv(k) r(k). (A2) 

Since the mass-to-light ratio of the accretion population is a 
constant factor and we are only interested in relative values, 
we can re- write Eq. Al as 



MAp(k)=fv(k)Lv(k), (A3) 

where both fv(k) and Ly(k) can be readily evaluated for 
each k shell. 

In order to obtain a realistic three-dimensional ra- 
dial profile, it is necessary to deproject the measured two- 
dimensional radial variation. The followed deprojection al- 
gorithm is based on the computation from outside to the 
center. In this sense, the first M\p(k) computed corresponds 
to the most external value (k = A^hdis, being A^hdis the to- 
tal number of concentric shells), which is not corrected for 
any projection effect. This initial value is then employed 
to obtain pAp(A r shciis) by dividing Map (Shells) by the vol- 
ume observed in projection in the external shell (in arbitrary 
units). In the successive steps we move towards the central 
regions (decreasing k), making use of the already determined 
PAp(i) values (with k < i < A r s h e iia). In addition 

^shells 

M A p(k)= J2 PAp(i)V(k,i), (A4) 

i=k 

where V(k, i) is the volume of the i th -shell in the line of sight 
of the fc th -shell, with k < i < A^heiis (see Fig. Al and figure 
caption) . Finally, we compute the desired spatial profile em- 
ploying 



/ Jv ahclls \ 

— Ul AP (k)~ £ PAp(i)V(k,i)j . (A5) 



Table 1. Observational configuration. 

Run 1 



Run 2 



Dates 
Telescope 
Spectrograph 
Detector 
Dispersion 
Wavelength Range 
Spectral Resolution 
Slit Width 
Spatial Scale 



9-10 Aug 1994 
CAHA 3.5m 

CTS 
CCD TEK 12 
3.46 A/pixel 
3700-7240 A 
8.1 A (FWHM) 

2.1 arcsec 
0.9 arcsec/pixel 



17-19 Dec 1995 

WHT 4.2m 
ISIS blue arm 
CCD TEK 1 
2.90 A/pixcl 
3740-6700 A 
12.3 A (FWHM) 

2.0 arcsec 
1.1 arcsec/pixel 
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Table 2. Galaxy sample. 
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Spectral gradients in central cluster galaxies 

Table 3. Central D4000 and Mg2 indices, and associated random errors (AD4000 and AMg2), for the central galaxies 
of the clusters named in the first column. The measurements correspond to a fixed metric aperture of 4 arcsec at the 
distance of the Coma cluster (~ 2.6 kpc). We have included the sample of CGA95, with revised central measurements 
for A 1795 and A 2199 (due to emission line contamination). The Mg2 values have been converted to the Lick system 
using the systematic offsets derived from the observation of template stars (see text). Mass deposition rate, M, and 
cooling radius, r coo i, for each cluster are taken from: [1] Peres et al. (1997, PSPC data), [2] Edge et al. (1992), [3] 
M092, [4] HBvBM, and [5] Arnaud (1988). When available, the values quoted for these last two parameters are 
the median, 10 and 90 percentile estimates. Central mass deposition rates (i.e. corresponding to the same metric 
aperture than that employed in the measurement of the line-strength indices) have been obtained assuming M oc r. 
The r coo i value for A 1126 is the mean cooling radii for the four galaxies with similar mass deposition rate (A 85, 
A 644, Hydra A and A 2597). 
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Table 4. D4000 and Mg2 gradients (dD4ooo/d log r, dMg 2 /dlogr), corresponding to the error-weighted least-squares fits shown 
in Figs. 3 and 4. Total gradients refer to fits computed by employing all the data excluding measurements within 1.5 arcsec 
and secondary nuclei. Inner and outer gradients correspond to fits inside and outside the emission line region, when present. 
Formal errors of the fits are given next to each gradient. 
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Spectral gradients in central cluster galaxies 17 

Table 5. Central D4000 indices derived from the nuclear (U— b) colours published by M092, using the 
linear relation fitted in Fig. 5b. Errors (AD4000) have been derived from the (U— b) errors reported 
by M092. Sources for mass deposition rate and cooling radius are those quoted in Table 3. 
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Figure 1: Comparison of our Mg2 measurements with 
those measured by the Lick group in a sample of 39 (run 1) 
and 13 (run 2) stars from the Lick library. The error bars 
in the lower right corner indicate the typical random er- 
ror. There is a systematic offset of 0.021 mag (dashed 
line) with a rms scatter of 0.010 mag. The offsets found 
analysing each run independently exhibit an excellent agree- 
ment: 0.020, rms = 0.008 and 0.021, rms = 0.021 mag for 
run 1 and 2, respectively. 

Figure 2: Nuclear spectrum of the central galaxy of 
the cluster Hydra A with typical emission lines. The arrows 
indicate the location of the Mg2 bandpasses. A simultane- 
ous fit of a scaled template spectrum, a 2 nd order polynomial 
and a Gaussian reproduces the local galaxy spectrum around 
the central bandpass. The final fit (polynomial+template) is 
shown as the thin line overplotted on the galaxy spectrum. 
The rms scatter of the residuals (shown in the inner small 
box) is 0.017, absolutely consistent with the mean random 
error in this spectral range (0.016). The pixels affected by 
the [Nl] A5199 emission are replaced by this combined fit. 
The [O in] A4959 emission in the blue bandpass is analo- 
gously removed. Common emission lines interpolated in the 
D4000 bandpasses are [Nem] A3869, [Sn] AA4069, 4076, and 
H<5. This procedure is repeated in each spectrum with emis- 
sion lines. 

Fig ure 3: D4000 gradients. See explanation in sec- 
tion [L2. 

Fig ure 4: Mg2 gradients. See explanation in sec- 
tion [l2. 

Figure 5: Panel (a): comparison of central D4000 mea- 
surements from JFN87 — open circles — and M089 — filled 
circles — with our central D4000 indices for four and three 
galaxies in common, respectively. Panel (b): correlation be- 
tween the nuclear (U— b) colours from M092 with the cen- 
tral D4000 values. 

Figure 6: Comparison of Mg2 gradients for four galax- 
ies in common with FFI95. Open symbols (squares: major 
axis, triangle: minor axis) are from FFI95, whereas filled 
circles correspond to measurements from CGA95 and this 
work (major axis, averaged at both sides of the galaxies). 
The minor axis data from FFI95 are rescaled to reproduce 
the same spatial scale than the major axis data. The up- 
per x-scale in the plots indicates the spatial scale in arcsecs. 
The atomic Mgb measurements of FFI95 were transformed 
into molecular Mg2 values and then converted to the Lick 
system by applying an offset of 0.014 mag. This offset has 
been derived from a comparison of the FFI95 data with a 
compilation of gradients by Gonzalez & Gorgas (1997). Our 
Mg2 indices were also transformed to the Lick system using 
the systematic offsets determined through the observation 
of template stars from the Lick library (see text). The short 
horizontal full lines at r ~ 1 arcsec in panels (a) , (b) and (d) 
correspond to the central Mg2 measurement given by Trager 
(1997). Similarly, the short horizontal dashed lines in pan- 
els (b) and (d) are the central Mg2 values measured by Lucey 
et al. (1997). The dotted arrows in the inner regions of the 
central galaxies of A 2199 and A 496 indicate the location 
of our Mg2 indices if we had measured Mgb in our spectra 
without removing the contamination by the emission line 
[Nl] A5199, and then transformed these Mgb values into 
Mg2 (in the same way as we did with the FFI95 data). The 
effect of the emission lines in A 2634 is completely negligible. 



Figure 7: Central D4000 and Mg2 measurements versus 
X-ray derived mass deposition rates (the sources for M are 
given in Table 3) . The central indices were obtained employ- 
ing a fixed metric aperture size of 4 arcsec projected at the 
distance of the Coma cluster. Central mass deposition rates 
have been computed from the total mass deposition rates as- 
suming M tx r. Error bars in mass deposition rate indicate 
the 10 and 90 percentile estimates. Filled and open sym- 
bols correspond to CCGs with and without emission lines 
in the central regions, respectively. Circles are galaxies from 
CGA95 and this paper, whereas triangles indicate D4000 val- 
ues derived from nuclear (U— b) colours published by M092, 
as explained in the text. 

Figure 8: Central D4000 versus Mg2 measurements. See 
explanation in section |5.l| . 

Figure 9: Central spectrum of three CCGs of our sam- 
ple with emission lines (A 1795, A 2597 and Hydra A, from 
top to bottom, -thick line-), compared with a scaled tem- 
plate spectra (-|-low-order polynomial) obtained from the 
central regions of CCGs without emission lines (A 2124, 
A 644 and A 644, from top to bottom, -thin line-). A 
marginal detection of He n A4686 is present in A 1795 and 
in A 2597, whereas this is uncertain in the case of Hydra A. 

Figure 10: Modelling of different star formation 
episodes overimposed on the spectrum of an old stellar 
population. We have employed the GISSEL96 predictions 
(Bruzual & Chariot 1997), with solar metallicity and Scalo 
(1986) IMF. Circles are relative UV fluxes from Crawford & 
Fabian (1993) —filled circles: A 1795, open circles: A 2597— 
, whereas the filled square is an IR flux estimation from the 
averaged nuclear (U-I) colours of A 1795 and A 2597 re- 
ported by M092. Panel (a) shows the spectral energy dis- 
tribution (SED) corresponding to an old (15 Gyr) stellar 
population (thick line), a continuous star formation during 
the last 10 Gyr (dotted line), and the final product of both 
stellar components (thin line) assuming that the fraction 
of V light that comes from the continuous star formation 
is fv = 0.3. The prediction of this continuous star forma- 
tion (thin line) is compared, in the following panels, with 
the final SEDs obtained by superimposing a single burst of 
star formation on the spectrum of the old stellar popula- 
tion. In particular, we have represented the resulting spec- 
trum (thick line = 15 Gyr stellar population + single burst, 
fv ~ 0.25) when 0.01, 0.1 and 0.2 Gyr have elapsed af- 
ter the time the burst took place (panels (b), (c) and (d) 
respectively). See discussion in the text. 

Figure 11: Comparison of the break radius obtained 
through the simultaneous fit of two straight lines forced to 
join at a common radius, with the radius of the emission- 
line region. Break radii and their error bars correspond to 
the mean and standard deviation derived in Monte Carlo 
simulations. When available, estimates of the break radius 
from D4000 and Mg2 line-strength gradients have been av- 
eraged. Emission-line radii where determined in the spec- 
troscopic images as the radii where the emission lines not 
longer detectable (see section 4.2). 

Figure 12: Line-strength gradients in D4000 and Mg2 
for the galaxies tabulated in Table 4 as a function of the 
mass deposition rate — panels (a) and (b). Open symbols 
are total gradients in galaxies without emission lines. Filled 
symbols correspond to inner gradients (those measured in 
the emission line region) . Error bars in mass deposition rate 
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indicate the 10 and 90 percentile estimates. In panels (c) 
and (d) we compare the mean gradients as a function of 
the galaxy type: normal ellipticals (Es), central dominant 
galaxies in clusters without cooling flow (CCGs no CF), and 
cooling flow galaxies without (CFGs no EL) and with emis- 
sion lines (CFGs EL). In the latter type we also distinguish 
between inner and outer gradients, which correspond to the 
emission line region and outside this region, respectively. 
The error bars in the last two panels are the rms scatter 
around the mean values. The mean D4000 gradient for Es 
(—0.315, rms = 0.090) has been computed from the data 
presented by Munn (1992) and Davidge & Clark (1994) (see 
CGA95 for further details). In addition, the Mg2 gradient 
for Es (—0.055, rms = 0.025) is the one derived by Gonzalez 
& Gorgas (1997) using published and unpublished data for 
109 galaxies with reliable Mg2 profiles. 

Figure 13: Total mass transformed into stars in the 
emission-line region of the central galaxies with reliable 
D4000 gradients in the outer parts. The plotted values have 
been derived using Eq. ^, and assuming that a single burst 
of star formation (after 0.1 Gyr) has produced the observed 
D4000 variation corresponding to the difference between the 
measured inner gradient and the extrapolated outer gradi- 
ent. 

Figure 14: Spatial profile of the mass transformed into 
stars obtained through the deprojection of the observed ra- 
dial variations in D4000 (see Appendix A). The dashed line 
is the expected density of the mass deposition profile is one 
assumes that M(< r) oc r. The numbers in parenthesis are 
the profile slopes obtained by least-squares linear fits. 

Figure Al: We have obtained the spatial profile of 
the density of mass transformed into stars pAp(fe) in each 
fc-shell by subtracting the contribution of the previously de- 
termined pap(*) (k < i < TVsheiis)- For this purpose it is nec- 
essary to evaluate the volume of each i-shell in the line of 
sight of the considered fc-shell, V(k,i). If we define v(r c ,r s ) 
as the volume of the intersection between a cylinder of ra- 
dius r c and a sphere of radius r s , it is easy to show that 
V(k, i) = v(r h ,r 2 ) - v(r a , r 2 ) - v(r h , n) + «(r a , n), where 
v(r c ,r s ) = 2nrlh+l-K{r s -h) 2 {2r s + h), being h 2 = r 2 - 
r 2 

1 c ■ 
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